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We investigate finite-size effects on tfie phase diagram of strong interactions, and discuss tfieir 
influence (and utility) on experimental signatures in high-energy heavy ion collisions. We calculate 
the modification of the pseudocritical transition line and isentropic trajectories, and discuss how 
this affects proposed signatures of the chiral critical endpoint. We argue that a finite-size scaling 
analysis may be crucial in the process of data analysis in the Beam Energy Scan program at RHIC 
and in future experiments at FAIR-GSI. We propose the use of extrapolations, full scaling plots 

. and a chi-squared method as tools for searching the critical endpoint of QCD and determining its 

' universality class. 
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■ I. INTRODUCTION 
^' 

, Although the process of phase conversion of hot hadronic matter that presumably happens in ultra-relativistic heavy- 
1 ion collision (HIC) experiments is often compared to the cosmological quark-hadron transition in the early universe, 

■ the relevant space-time scales differ by almost twenty orders of magnitude. In particular, while one would need 
temperatures T ~ 10^ GeV to be subject to appreciable finite-size effects in the case of the early universe, which 
of course spares the primordial quark-hadron transition from these effects, this is not the case for the quark-gluon 
plasma formed in HICs [J, Therefore, in the context of HICs, descriptions of an equilibrium quark-gluon plasma 
in the thermodynamic limit are, in principle, missing non-negligible finite-size corrections. 

^ I Strictly speaking, from the statistical mechanics standpoint, there is no phase transition in a system of finite 
size, and no criticality can be defined. Singularities that would play the role of strong lighthouses in the search for 
criticality are rounded and become smooth. Actually, the originally critical quantities are also shifted in the T — fi 
plane according to volume variations, becoming pseudocritical, with some unavoidable ambiguity in their definition. 
This is the actual scenario in any experiment in HICs. In addition, the time scale involved in the collision process is 
also not large compared to the microscopic time scale to achieve equilibrium. Therefore, thermodynamic quantities 
directly inferred from data assuming an infinite system in equilibrium should not correspond to its true and unique 
\ value in the thermodynamic limit. Moreover, even neglecting dynamical non-equilibrium effects, all signatures of 
the second-order critical endpoint based on non-monotonic behaviour 0, 0] (or sign modifications Q) of particle 
correlation fluctuations will probe a pseudocritical endpoint, shifted significantly from the genuine critical endpoint 
by finite-size corrections and boundary effects. 
' In this paper, we do not address non-equilibrium effects, but argue that finite-size modifications will play a crucial 
, role in the search for the critical endpoint in HICs. On one hand, they restrict the access of the proposed particle 
correlation fluctuation signatures to (volume- and boundary-condition-dependent) pseudocritical quantities. On the 
other hand, the fact that the experimental data can be arranged as a set of systems of different sizes can actually help 
us extract information directly related to the thermodynamic limit of QCD, especially concerning its critical endpoint 
[ and the universality class of the transition. 

First, we investigate finite-size effects on the phase diagram of strong interactions, presenting results for the mod- 
ification of the pseudocritical line in the temperature-chemical potential phase diagram, as well as for isentropic 
trajectories. For this purpose, we employ the linear sigma model coupled to quarks with two flavors, Nf — 2 6]. This 
effective theory has been widely used to describe different aspects of the chiral transition, such as thermodynamic 
properties 0, S S [3 ^^'^ the nonequilibrium phase conversion process ^llj , as well as combined to other models in 
order to include effects from confinement We show that the volume-dependence of the QCD phase diagram in 
the regime of energy scales probed by current heavy ion experiments can be large. 

Second, treating HIC data separated into slices of different impact parameter width as a realization of a set of 
thcrmodynamical systems of different sizes and temperatures, we propose the use of full scaling plots of cumulants of 
fluctuations of pions as a tool for searching the critical endpoint of QCD and determining its universality class. From 
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this point of view, HICs could be seen as the "experimental realization" of lattice simulations for systems of different 
sizes. In this picture, signatures of the second-order critical endpoint based on non-monotonic behaviour (or sign 
modifications) of particle correlation fluctuations will actually probe a set of pseudocritical endpoints, a feature which 
could contribute to the broadening of the experimental signal, helping to wash it out in the large thermal background. 
As in Monte Carlo simulations, the peculiarities of the finite system become relevant (technically, this corresponds to 
the need of inclusion of irrelevant operators in the formalism). In that case, one has either to use the knowledge of 
the peculiarities of the system or resort to finite-size scaling (FSS) techniques to extract indirect information about 
criticality. 

Finally, we also propose a chi-squared method for searching the critical endpoint of QCD and determining its 
universality class that could be useful in the process of data analysis in HIC experiments, in particular in the Beam 
Energy Scan program at RHIC-BNL and at FAIR-GSI. 



II. FINITE-SIZE CORRECTIONS TO THE PHASE DIAGRAM OF QCD 

Finite-size effects in the study of phase transitions via lattice simulations, seen as an inevitable drawback of the 
method, and the necessary extrapolation to the thermodynamic limit have been thoroughly studied for decades p^ . 
the ultimate answer to the problem being given by the method of FSS [T3] . Systematic calculations of finite- volume 
corrections can then be computed not only near the criticality of continuous (second-order) transitions, but also for 
first-order phase transitions [l5l |. so that the thermodynamic limit can be taken in the calculation of properties of 
phase diagram. 

However, for natural systems that are truly small, one should study the modifications caused by its finitude in the 
phase diagram before comparing to experimental observables. This is the situation in the study of the chiral transition 
in HICs. Contrasting to the large number of studies in the computation of the chiral condensate and related quantities 
within chiral perturbation theory (see e.g. [3 and references therein), this issue is often overlooked in the case of the 
quark-gluon plasma. Amon g th e exceptions, there is a lattice estimate of finite-size effects in the process of formation 
of the quark-gluon plasma [l7| , a few studies within the Nambu-Jona-Lasinio model , and an analysis of color 
superconductivity in dense QCD in a finite volume [l^, besides investigations of finite-size effects on the dvnamicsfsij 
of the plasma [U, |2l . The effect of finite volume on the pion mass and chiral symmetry restoration at finite temperature 
and zero density was investigated by Braun et al. using the proper-time renormalization group [20j . 

To describe the chiral phase structure of strong interactions at finite temperature and baryon density, we adopt the 
linear sigma model 

C = [i^d^ + /i7° - 9^)] + \di.od^o - y{o, 7?) , (1) 
where /i is the quark chemical potential, 

V{a,^)^-^{a'^-v''f-ha (2) 

is the self- interaction potential for the mesons, exhibiting both spontaneous and explicit breaking of chiral symmetry. 
We dropped the pion terms for simplicity [s^]. The Nj = 2 massive fermion fields represent the up and down 
constituent-quark fields V' = {u,d). The scalar field a plays the role of an approximate order parameter for the chiral 
transition, being an exact order parameter for massless quarks and pions. Quarks constitute a thermalized fiuid that 
provides a background in which the long wavelength modes of the chiral condensate evolve. We incorporate quark 
thermal fluctuations in the effective potential for cr, i.e. we integrate over quarks to one loop. The parameters of the 
lagrangian are chosen such that the effective model reproduces correctly the phenomenology of QCD at low energies 
and in the vacuum. So, we impose that the chiral 51^^(2) (g) S'[/fl(2) symmetry is spontaneously broken in the vacuum, 
and the expectation values of the condensates are given by (cr) — and (tt) = 0, where = 93 MeV is the pion 
decay constant. The explicit symmetry breaking term is determined by the PCAC relation which gives h = /jrmj, 
where ~ 138 MeV is the pion mass. This yields — f"^ ~ rn^/X. The value of A = 20 leads to a cr-mass, 
TO^ = 2A/^ -|- toJ, equal to 600 MeV. The choice g = 3.3, which yields a reasonable mass for the constituent quarks 
in the vacuum, Mq = gf^, leads to a crossover at finite temperature, T, and vanishing chemical potential, so that 
the first-order transition line that starts at vanishing temperature and nonzero chemical potential stops at a critical 
endpoint. 

In the case of a finite system of linear size L, the momentum integral from the one-loop quark contribution to the 
effective potential is substituted by a sum 

^ = ^ E h (l + ^-(^-^^/^) + log (l + , (3) 
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where E]^ = i/k^ 



and 



— g\(j\ is the effective mass of the quarks 33|. It is clear that finite-size effects 



in a quantum field theory are, in a way, very similar to thermal effects with 1/L playing the role of the temperature 
[2l| . However, in the latter boundary conditions are determined by the spin-statistics theorem, whereas in the former 
there is no clear guidance (see, e.g., discussions in Refs. [H, [23] )■ Choosing periodic boundary conditions (PBC) 
are useful to focus on bulk properties of the model undisturbed by surface effects (111, , but one is always free 

to choose anti-periodic boundary conditions (APC) or any other, instead. Depending on the physical situation of 
interest, a boundary condition may be more realistic than others. To illustrate the dependence of the pseudocritical 
phase diagram on the boundary conditions, we show results for PBC and APC js^]. For PBC, e.g., the components 
of the momentum assume the discretized values = 2'Kii/L, being integers, and there is a zero mode (absent 
for APC, for which ki = Tr{2£i + 1)/L). At zero temperature, this zero mode in the case with PBC will modify the 
classical potential, generating size-dependent effective couplings and masses. 



APC 




L = infinity 

- ■ L = 5 fm 

L = 3 f m 

■ ■ L = 2fin 



200 



300 400 500 

H ( MeV) 



600 




300 400 

\l ( MeV) 



FIG. 1: Pseudocritical transition lines and endpoints for different system sizes within the linear a model with periodic (left) 
and antiperiodic (right) boundary conditions. 



In what follows, we analyze systems of sizes between 10 fm and 2 fm which relate to the typical linear dimensions 
involved in central and most peripheral collisions of Au or Ph ions at RHIC and LHC, respectively. Figure [T] displays 
the shift of the pseudocritical transition lines and their respective endpoints as the size of the system is decreased. 
The transition lines represent pseudo- first-order transitions, characterized by a discontinuity in the approximate order 
parameter, the chiral condensate cr, and the production of latent heat through the process of phase conversion. We 
find that those lines are displaced to the region of higher /x and shrinked by finite-size corrections. The former effect 
is sensibly larger when PBC are considered, indicating that the presence of the spatial zero mode tends to shift the 
transition region to the regime of larger chemical potentials. Both boundary conditions reproduce the infinite- volume 
limit for L > 10 fm. Figure [5] shows the corresponding displacement of the pseudocritical endpoint, comparing PBC 
and APC: both coordinates of the critical point are significantly modified, and /icep is about 30% larger for PBC. 
For fjL = Q, the crossover transition is also affected by finite-size corrections, increasing as the system decreases, as 
shown in Figure [Til Again, PBC generate larger effects: up to ~ 80% increase in the crossover transition temperature 
at ^ = when L = 2 fm. 




FIG. 2: Displacement of the pseudocritical endpoint in the FIG. 3; Normalized crossover temperature at /i = as a 
r — /i plane as the system size is decreased for different function of the inverse size 1/L for the cases with PBC 
boundary conditions. and APC. 
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Results for the isentropic trajectories are shown in Figure [H comparing the infinite- volume limit with the finite 
system with L = 2 fm in the cases with PBC and APC. For sufficiently high temperatures, the isentropic lines in 
the thermodynamic limit are reproduced, while large discrepancies are found around and below the transition region. 
The remarkable variations at low temperatures are due to the shift of the pseudocritical line and, for PBC, to the 
zero-mode-induced modifications of the vacuum properties, especially the vacuum constituent quark mass. Although 
there is no strong focusing effect around the critical endpoint (as observed previously for this chiral model Q and 
similar ones [26l]), it is clear that the density of isentropic trajectories around the pseudocritical endpoint is sensibly 
higher than in the thermodynamic limit. 
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FIG. 4: Isentropic trajectories, labeled by the respective value of entropy per baryon number, in the thermodynamic limit 
(solid, red lines) and for L = 2 fm with PBC (left) and APC (right). The red dot is the genuine critical endpoint, while the 
square is pseudocritical one. 



III. SCALING PLOTS AND THE UNIVERSALITY CLASS OF QCD 

The close neighborhood of critical points exhibit strong and peculiar scaling behavior that can be used as signals of 
its presence [l3, 14] . This fact has been recognized long ago as a possible means of detecting experimentally the critical 
endpoint of QCD by using event-by-event analysis in HIC experiments Q. In fact, higher moments of fluctuations of 
particle multiplicities diverge with increasing powers of the correlation length Q . As discussed in the introduction and 
shown quantitatively for the chiral phase diagram in the last section, finite-size corrections are relevant in the context 
of HICs, not only smoothcning those singularities but also generating shifted peaks and pseudocritical observables. 
Based on the analysis of the last section, the amount of displacement of the pseudocritical transition parameters that 
are probed at each experimental event should present a significant centrality (size) dependence. Therefore, when 
considering data from a centrality window, there is an average between different peaks which should broaden even 
further the non-monotonic signal of the pseudocritical endpoint. Furthermore, if the transition line is indeed strongly 
shrinked and shifted to higher chemical potentials, the current experiments may not be able to probe this regime of 
pseudocritical peaks and no non-monotonic behaviour would be found in the data. Finally, if this is not the case, to 
obtain, within this method, experimental information about the genuine QCD critical endpoint, one would need to 
gather a set of pseudocritical endpoint measurements for different system sizes and apply a procedure of extrapolation 
to the thermodynamic limit, namely a FSS analysis. 

Recently, it was proposed that a FSS analysis on top of the critical point would be a more efficient way of searching 
[2?! , and, as mentioned previously, FSS was also used before in the context of the dynamics of a first-order transi- 
tion in QCD in a finite- volume quark-gluon plasma [l|, Q • The caveat of the aforementioned method for searching 
experimentally the critical endpoint of QCD is the assumption that one is essentially on top of the criticality. This 
condition will be hardly met in HIC experiments, at least with enough statistics to stand out of the background. Even 
if one is very close, or on top, of the critical point, a system that is as small as the quark-gluon plasma created under 
these conditions, where all relevant length scales are close 3, will be strongly affected by finite-size effects. Namely, 
all singularities will be rounded and spread, something that can provide misleading artifact "signals" even in a FSS 
analysis when restricted to the assumption of being on top of the critical point. 

The alternative is simply making use of the full power of FSS, which is a method that is valid also when one is 
not so close to the critical point [35|, provided that we use full scaling plots. This technique is predictive even for tiny 
systems in statistical mechanics (see, e.g., Ref. [1^), and, since it can be used in various different regions of the phase 
diagram not far from the critical point, it can in principle provide enough statistics for data analysis. We propose the 
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use of this method to search for the critical point of QCD and determine its universahty class, making profit out of 
the fact that we can build out of heavy ion data a set of systems of different sizes. 

The FSS hypothesis [3, 0, [2§| , that can be derived applying the renormalization group to critical phenomena 
[U, [3^ , states that an observable X in a finite system at temperature T can be written, in the neighborhood of 
criticality, in the form [24| 

X{t,g-i-L)^L'^^/'^f,{tL^/'^) , (4) 

where t = (T — Tc)/Tc is the reduced temperature (Tc being the temperature associated with the critical point), 
is the bulk (dimension) exponent oi X, g the dimensionless coupling constant, I a length to fix the renormalization 
procedure, and v the critical exponent of the correlation length that diverges as ^ '--^ at criticality 36] . The function 
fx{y) is universal up to scales fixing, and the critical exponents are sensitive essentially to dimensionality and internal 
symmetry, which will give rise to the different universality classes 1131 • 

The critical contribution to experimental observables such as moments of fluctuations of particle multiplicities near 
the critical point ^, i^J will suffer sizable corrections from the function fx{y 7^ 0) that can be computed within a 
model. One should note that the observables satisfying the scaling relation ( [4]) are the ones directly related to 
correlation function of the order parameter of the transition, such as fluctuations of the multiplicity of soft pions, as 
discussed in Ref. Therefore, a careful investigation of how the evolution of the system in the hadronic phase and 
the large background of thermal correlations might affect or even hide the scaling behaviour of a given experimental 
observable is demanded. Looking for the collapse of points representing different sizes of the system at < = [23| can 
be dangerously misleading when there is not enough statistics as in HICsjSTj. 

Plotting observables such as cumulants of fluctuations of soft pions properly normalized [s^ and as a function of 
the full scaling variable tL^/'^ can be an efficient method for searching the critical endpoint of QCD and determining 
its universality class. The critical temperature and chemical potential (the coordinates of the critical endpoint) as 
well as the critical exponents would be parameters in a fit (scaling plot), in the same manner as is done in statistical 
mechanics simulations. In HICs, one can arrange the data as representing systems of different sizes by using centrality 
bins, or by directly using the number of participants in the collision, instead of L, in the scaling plot[3^. or each 
system of a given size, and also a given center-of-mass energy -y/s, one can estimate the temperature and chemical 
potential employing thermal models in the chemical freeze-out, or use directly 5s = {^fs — ^/s^) / ^/s^, y/s^ substituting 
(Tcfic) as the free parameter associated with the position of the genuine critical endpoint. When one plots XL^^^^" 
as a function of, e.g., T (not a scaling plot), one finds scattered curves that cross in the critical point. In a full scaling 
plot, though, that uses the correct scaling variable, all these curves collapse, provided that the system is not far from 
the critical point, as a consequence of scaling and universality. By finding the correct critical exponents in the case 
of the QCD phase transition in HICs, one would determine its universality class. 

One can search for this full scaling behavior in different cumulants of fluctuations and in different parts of the phase 
diagram, increasing dramatically the statistics. Still, realistically, the thermal background of particle production could, 
in principle, hide the critical contribution inside the error bars and the process of flnding the correct value of the 
critical exponent could be difficult. Therefore, to put in practice the idea of scaling-plot fits described above, we may 
introduce the following chi-squared method. Once an ensemble of experimental results for the properly normalized 
observable is obtained, one can find the best scaling fit via the minimization of the quadratic deviation among them. 
For example, let 

Q (Vs. ^)=X N;2t , (5) 

for the observables X (say, cummulants of soft pion fluctations). Then we may define the following function 

X {v, V^; y) = 2^ -— - 1 

^part ■5s = a 

where the sums are taken over all pairs of (\/s, -y/A^part) keeping the scaling variable y constant, with 

If all of X values for different (^/s, i/iVpart) obey the finite size scaling, then the above should vanish at the correct 
value of V and ^/s'c for any y. For energies to which the observables attain values close to the critical point, then we 
should have a very week y-dependence of v and ^fsl and also the value of at the minimum, only coming from the 
statistical errors. In practice, we plot for a given y (or various combinations of (\/so' ■\/-^part,o) and localize the 
minimum in v and 
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Thus, by studying the energy dependence of the minimum of the above quantity, we may identify the energy value 
which passes close to the critical point with the free variational parameters being the critical exponent v and the 
critical center-of-mass energy ^fsZ (and 7^ in the general case). 

IV. FINAL REMARKS 

In this paper, we have discussed different aspects of finite-size effects that render crucial features in various proposed 
signatures of the QCD critical endpoint in HICs. On one hand, in addition to the usual smoothening and broadening 
of singularities, the finiteness of the physical system under a phase conversion process shifts pseudocritical observables 
associated with pseudocritical parameters with respect to the genuine criticality, in the thermodynamic limit. We have 
shown that those displacements in the case of the chiral transition can be large for energy scales typically encountered 
in current experiments, shrinking the T-/i domain where discontinuities of the chiral condensate are found and pushing 
the pseudocritical lines further into the high chemical potential region. These findings should affect all signatures of 
the critical endpoint based on non-monotonic behaviours or sign changes in particle correlation fiuctuations, probably 
contributing to the broadening of the signal within a not sufficiently small centrality window. 

Based on this considerations we proposed the use of the full power of the FSS hypothesis as a tool to detect and 
locate the genuine chiral critical endpoint of QCD, and determine its corresponding universality class, using HICs 
experiments. We also described a chi-squared method that can be applied to data in a systematic manner to investigate 
the presence of the scaling phenomenom in the Beam Energy Scan project at RHIC and in future observations at 
FAIR. 

Finally, to address the question of the viability of these different proposed signatures of the QCD critical endpoint 
when immersed in a large hadronic thermal background, as the one expected in the case of HICs, an interesting test 
to make would be the comparison of Monte Carlo simulations of thermal ensembles with and without a small fraction 
of critically correlated soft pions. 
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